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ABSTRACT 



We use XMM-Newton data to carry out a detailed study of the Si, Fe and Ni abundances in the cool cores of a representative sample 
of 26 local clusters. We have performed a careful evaluation of the systematic uncertainties related to the instruments, the plasma 
codes and the spectral modeling finding that the major source of uncertainty is in the plasma codes. Our Si, Fe, Ni, Si/Fe and Ni/Fe 
distributions feature only moderate spreads (from 20% to 30%) around their mean values strongly suggesting similar enrichment 
processes at work in all our cluster cores. Our sample averaged Si/Fe ratio is comparable to those measured in samples of groups 
and high luminosity ellipticals implying that the enrichment process in ellipticals, dominant galaxies in groups and BCGs in clusters 
is quite similar. Although our Si/Fe and Ni/Fe abundance ratios are fairly well constrained, the large uncertainties in the supernovae 
yields prevent us from making a firm assessment of the relative contribution of type la and core-collapsed supernovae to the enrichment 
process. All that can really be said with some certainty is that both contribute to the enrichment of cluster cores. 
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1. Introduction 



The intra-cluster medium (ICM), that fills the deep potential well 
of galaxy clusters, is rich in metals. In the high mass regime of 
rich clusters (M VJ > > 1O I4 M ), where temperatures and densi- 
ties of the gas are of the order of 3 - 10 keV and 10~ 3 cm -3 , 
respectively, the ICM is heavily ionized and the main processes 
producing emission lines are collisions of ions with electrons 
in a plasma in collisional ionization equilibrium. Since the el- 
ements are well confined within the cluster potential well they 
accumulate over the whole cluster history and retain important 
information on cluster formation and evolution. For gas temper- 
atures > 3 keV the most prominent emission lines are from He- 
like and H-like K-shell transitions of Fe (at ~ 7 ke V), along with 
less pronounced K-shell lines of S, Si (at ~ 2 keV), Ar, Ca and 
Ni (at ~ 8 keV). Temperature and abundances of the ICM are 
measured at the same time from the X-ray spectrum: the temper- 
ature is measured from the continuum emission that is almost en- 
tirely given by thermal bremsstrahlung, whereas the abundance 
of an element is derived by measuring the equivalent width of the 
line (once the continuum is known), that is directl y proportional 
to the ion-to-Hydrogen concentrations ratio (see Kaastra et al. 
2008, for a review of the thermal radiation processes in clus- 
ters). While in principle the abundance measurement is straight- 
forward, in practice various sources of uncertainties are present, 
amongst which (1) the accuracy of the atomic physics, (2) the 
moderate spectral resolution of the current imaging instruments 
which often results in line blending, this is particularly severe 
for L-shell blends where the transitions are not entirely under- 
stood, (3) the presence of temperature gradients in the ICM, es- 
pecially in the cluster cores, that need specific spectral modeling. 
With the advent of Chandra and XMM-Newton satellites, carry- 
ing detectors with both high spatial and spectral capabilities, the 
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statistical quality of cluster spectra (in particular for cool-core 
regions) improved dramatically. While on the one hand the sta- 
tistical errors associated with the abundance measurements have 
greatly decreased, on the other surprisingly little attention has 
been devoted to systematic sources of uncertainties which, un- 
der these circumstances, are likely to play an important role. 

In this paper our first goal is to provide a robust measure- 
ment of the distribution of abundances and abundance ratios of 
the chemical elements in the cores of a representative sample of 
nearby and bright cool-core clusters. In this context "robust esti- 
mate" essentially means that we will include in the error budget 
a careful evaluation of the systematic uncertainties potentially 
affecting our data. We will consider both instrument related sys- 
tematics and plasma code systematics. We stress that in this pa- 
per we are interested in the measure of the global abundances 
in the central region of each cluster, not in the analysis of radial 
profiles. Moreover these central regions will provide us with the 
maximum photon statistics because of their very intense central 
surface brightness peaks, thereby allowing us to explore the most 
subtle sources of systematic errors and reliably measure the most 
abundant elements observable in the ICM (e.g. Si, Fe and Ni). 

The elements observed in the ICM are produced through 
thermo-nuclear nucleo-synthesis in superno vae (SN) explosions 
occurring in the mem ber early-type galaxies (lArnaud et al.l l992; 
iRenzini et al.l 119931) . and they are eventually ejected/diffused 
into the ICM through ga lactic winds (for a review see 
Bland-Hawthorn et al. 2007) and ram pressure stripping (e.g. 
ISchindler & Diaferioll2008l and references therein). Supernovae 
are classified on the basis of their progenitor models. Nowaday 
two main groups of SNe are known: type la supernovae (SNIa) 
that derive from an accreting white dwarf in a binary system, 
and, core-collapsed supernovae (SNcc) whose progenitors are 
single massive stars (> 8M Q ). One of the most important quan- 



2 



S. De Grandi & S. Molendi: Metal Abundances in the Cool-Cores of Galaxy Clusters 



tity that theoretical models of SNe provide are the yields of 
the various chemical elements, namely the mass per element 
and per SN event. While, SNIa produce mostly Fe, Ni and Si- 
group metals (i.e. Si, S, Ar, Ca), the SNcc ejecta are rich in a- 
elements (i.e. O, Ne, Mg, as well as few Si-group elements). 
In the 90's BeppoSAX and ASCA observations revealed that the 
centers of cool-core clusters always display an excess (with re- 
spect to the extern al cluster regions) of Fe around the central 
cluster galaxy (e.g.[F ukazawa et alj2000t iDe Grandi & Molendil 
2001; Leccardi & Molendi 2008a). The iron mass inferred from 
this excess is consisten t with being entirely produced by the 
giant galaxy itself (e.g. IDe Grandi et al.l l2004t iBohringer et alj 
2004), which is always found in these systems. Other ele- 
ments suc h as Mg, Si and S show abundance peaks in cluster 



cores (e.g. Fukazawa et al. 19981 : Finoguen ovet al.ll2000l [2001; 
Tamu ra et alJ l2004t ISato et alJ l2009a)~. Under the assumption 
that the sole source of metals are the two types of SNe, ob- 
servations of the a-elements/Fe abundance ratios coupled with 
the knowledge of the SNe chemical yields allow the determi- 
nation of the SNIa and SNcc proportion within each cluster. 
Finoguen ov et alJ (12000) found an increasing Si/Fe ratio with ra- 
dius in clusters indicating a greater predominance of SNcc en- 
richment at large radii, while the innermost parts appeared dom- 
inated by SNIa products. However works based on data taken 
with the more recent XMM-Newton and SUZAKU observato- 
ries have not fully confir med this overall picture for rich clus- 
ters. [Tamura et al. (2004) studying the spatial distributions of 
metals in a sample of cool-core clusters observed with XMM- 
Newton found uniform Si/Fe and S/Fe ratios within the ICM, 
but increasing O/Fe (although this ratio was prone to large un- 
certainties). SUZAKU observations confirm these finding s (e.g. 
ISato et alj|2007al: iMatsushita & Suzaku SWG teamll2008l) . The 
relative proportion of SNIa and SNcc found by XM M-Newton 
and S UZAKU observations are in raw agreement (Sat o et all 
2007a), although these results were achie ved under the choice 
of specific compilat ions of SNe yields (e.g. Iwa moto et al. 1999; 
iNomoto et al.ll2006h . It is worth noting that none of the combi- 
nations of the theoretical SNIa and SNcc yields available in the 
literature have reproduced the overall elemental pattern in clus- 
ter core s (fits were not formally acceptable based on the y 2 val- 
ues; e g. lBaumgartner etaT]|2005t Ide Plaa et ail2007l ISato et alJ 
2007a). There are at least two possible solutions to this prob- 
lem: the existence of an addit ional source of elements other than 
SNIa and SNcc (seelFinoguenov et al.ll2.Q02l: iBaumgartner et alJ 
120051: iMannucci et al.l2006l). or. the need of a rev ision of the the- 
oretical SNe models (e.g. Young & Fryer 2007). Interestingly a 
systematic exploration of the full rang e of the uncertaintie s of 
theoretical yields was last provided by iGibson et al.l d!997l) . al- 
though since then, sever al new compilations of SNcc yields have 
becom e available (e.g. IChieffi & Limongil 120041: INomoto et alJ 
2006) along wi th new yields for SN Ia dlwamoto et alj|1999l) . As 
pointed out by IGibson et al.l d 19971) the relative contribution of 
SNIa and SNcc to the Fe abundance in the ICM depends cru- 
cially upon the adopted theoretical SNe models. 

The second main goal of this paper is to provide a critical as- 
sessment of the relative contributions of the SNIa versus SNcc to 
the enrichment process. C ontrary to other sim ilar works present 
in the recent literature (e.g. Ide Plaa et al.l2007l) we do not neglect 
the uncertainties associated to the current theoretical yields and 
explore how they affect the derived SNIa fraction. 

This paper is organized as follows. In Sect. 2 we present the 
sample, in Sect. 3 we describe the cleaning process of the raw 
XMM-Newton archival data, and, in Sect. 4, we concentrate on 
the spectral analysis and on the choice of the best fitting model. 



In Sect. 4 we also discuss the systematic uncertainties related 
to imperfections in the cross-calibration between the three EPIC 
instruments. In Sect. 5 we present results on the derived distri- 
butions of metal abundances and abundances ratios. In this same 
section we compare results from two different choices for the 
spectral codes. In Sect. 6 we discuss the relative contribution to 
the overall enrichment process of different SNe types. In Sect. 7 
we discuss our findings and in Sect. 8 we summarize our main 
results. 

We assume Hq = 70 km s Mpc" 1 and Qa = 0.7, the source 
redshifts are all extracted from the NASA/IPAC Extragalactic 
Database (NED). All metallicity measurements we show in 
this paper are relative to t he photospheric solar abundances of 
lAnders & Grevessel dl989l) . In this set of abundances, Fe has a 
number density of 4.68 x 10" 5 , Si 3.55 x 10 s and Ni 1 .78 x 10~ 6 , 
all number densities are relative to H. All uncertainties shown 
are l<x confidence level (i.e 68%). 



2. The sample 

St arting from t he flu x-limited sample of 55 galaxy clusters listed 
bv lEdge et all (1 19901) and the ROSA T PSPC analysis of this sam- 
ple performed by |Peres"e t al. ( 1998), we have selected all the ob- 
jects r wjfh_amass deposition rate different from zero (see Table 
5 in iPeres et al.lll998l) . Among these clusters we have chosen 
all those that were observed with XMM-Newton EPIC within 
April 2008, with the exclusion of clusters whose observations 
are strongly affected by soft protons (e.g. Abell 644 and Abell 
2142). The 32 clusters that meet these requirements are listed in 
Table Q] 



3. Data preparation 

We reprocessed the Observation Data Files (ODF) using the 
Science Analysis System (SAS) version 7.0.0. After the produc- 
tion of the calibrated event lists for the EPIC MOS1, MOS2 and 
pn observations with emchain and epchain tasks, we have per- 
formed a soft proton cleaning using a double filtering process. 

We first removed soft protons spikes by screening the light 
curves produced in 100 bins in the 10-12 keV band and applying 
an opportune threshold for each instrument, generally a thresh- 
old of 0.20 cts s- 1 forMOSl and MOS2, and of 0.60 cts s 1 for 
pn. Then to eliminate possible residual flares contributing be- 
low 10 keV, we extracted a light curve in the 2-5 keV band and 
fitted the histogram produced from this curve with a Gaussian 
distribution. To generate the final filtered event file we have re- 
jected all events registered at times when the count rate was 
more than 3cr from the mean of this distribution. Finally, we 
have filtered events files according to FLAG (FLAG==0) and 
PATTERN (PATTERNS 12 for MOS and PATTERS==0 for pn) 
criteria. The resulting effective exposure times of the observa- 
tions are reported in Table Q] 

Using the cleaned events files we have extracted spectra from 
a circular region centered on the clusters emission peaks. The 
physical radius of this r egion is 0.5 r cooi, w here r coo i is the cool- 
ing radius computed bv lPeres et al.l (Il998l) . With this choice of 
the extraction radius we sample a significative portion of the cool 
core and we assure that the extracted region is contained within 
the EPIC field of view for all our clusters. Prominent point-like 
sources ha ve been removed from the extraction region. Othe r au- 
thors (e.g. Ide Plaa e t al. 2007; iRasmussen & Ponmanll2007l) se- 
lect there extra ction regions as fixed fractions of the cluster scal- 
ing radius, e.g. Ide Plaa et al.1 (120071) used 0.2 of rsoo, where rsoo 
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Table 1. Starting sample of 32 galaxy clusters extracted from the X-ray flux limite d B55 sample (Ed ge et al]|1990h . All objects are 
cool-core clusters with a mass deposition rate different from zero dPeresetalJI 19981) . 



Cluster 


ObsID 


t 


"h 




Exp.time c 




''coo//2 d 








10 2 W 2 


MOS1 


MOS2 


pn 


arcmin (kpc) 


Abell 85 


0065140101 


0.0551 


0.0342 


12.5 


12.5 


9.7 


1.137 (73) 


Abell 262 


0109980101 


0.0163 


0.0538 


23.8 


23.9 


19.2 


2.613 (52) 


AWM7 


0135950301 


0.0172 


0.0983 


30.7 


30.7 


28.7 


2.473 (52) 


Abell 3112 


0105660101 


0.0750 


0.0261 


23.3 


23.3 


18.6 


1.124(96) 


Perseus 


0085110101 


0.0179 


0.1410 


53.7 


53.7 


51.3 


4.263 (93) 


2A 0335 + 096 


0147800201 


0.0349 


0.1780 


80.2 


80.9 


73.9 


2.596 (108) 


Abell 478 


0109880101 


0.0881 


0.1510 


28.7 


79.0 


42.6 


1.032 (102) 


Abell 496 


0135120201 


0.0329 


0.0458 


17.8 


17.6 


14.1 


1.397 (55) 


PKS 0745 - 191 


0105870101 


0.1028 


0.4240 


18.9 


18.9 


9.8 


0.944 (107) 


Hydra A 


0109980301 


0.0549 


0.0494 


18.6 


19.8 


15.8 


1.266 (81) 


Abell 1060 


0206230101 


0.0126 


0.0490 


40.0 


41.5 


33.2 


2.588 (40) 


Virgo (M87) 


0200920101 


0.0044 


0.0254 


77.8 


79.4 


70.0 


9.352 (51) 


Centaurus 


0046340101 


0.0114 


0.0806 


40.7 


41.3 


40.3 


2.928 (41) 


Abell 1644 


0010420201 


0.0473 


0.0499 


15.1 


15.2 


12.8 


0.521 (29) 


Abell 1650 


0093200101 


0.0838 


0.0156 


38.1 


36.6 


34.4 


0.878 (83) 


Abell 1651 


0203020101 


0.0849 


0.0181 


10.2 


10.8 


6.9 


0.669 (64) 


Abell 1689 


0093030101 


0.1830 


0.0182 


37.5 


37.6 


32.6 


0.520 (96) 


Abell 3558 


0107260101 


0.0480 


0.0389 


44.1 


44.0 


38.9 


0.602 (34) 


Abell 3562 


0105261301 


0.0490 


0.0385 


40.3 


40.2 


38.0 


0.834 (48) 


Abell 3571 


0086950201 


0.0391 


0.0371 


24.3 


24.5 


16.5 


1.119(52) 


Abell 1795 


0097820101 


0.0625 


0.0119 


39.5 


39.6 


22.7 


1.232 (89) 


Abell 2029 


0111270201 


0.0773 


0.0304 


10.7 


11.2 


10.3 


1.059 (93) 


Abell 2052 


0109920101 


0.0355 


0.0272 


30.3 


30.3 


26.5 


1.747 (74) 


MKW3s 


0109930101 


0.0450 


0.0303 


38.7 


38.3 


34.1 


1.620 (86) 


Abell 2065 


0202080201 


0.0726 


0.0295 


20.7 


20.7 


16.6 


0.338 (28) 


Abell 2063 


0200120401 


0.0349 


0.0298 


8.1 


8.9 


5.4 


1.152(48) 


Abell 2199 


0008030201 


0.0302 


0.0086 


14.6 


14.6 


10.0 


1.985 (72) 


Abell 2204 


0112230301 


0.1523 


0.0567 


19.9 


20.2 


13.7 


0.629 (100) 


Cygnus A 


0302800101 


0.0561 


0.3490 


22.2 


22.2 


20.2 


1.041 (68) 


Abell 2597 


0147330101 


0.0852 


0.0249 


52.2 


53.3 


48.5 


0.792 (76) 


Abell 4038 


0204460101 


0.0300 


0.0155 


29.4 


29.2 


27.6 


1.859 (67) 


Abell 4059" 


0109950101 


0.0475 


0.0110 


12.7 


13.7 


7.0 


1.378 (77) 


Abell 4059 


0109950201 


0.0475 


0.0110 


23.3 


23.2 


19.7 


1.378 (77) 



Notes: 5 redshifts taken from the NASA Extragalactic Database; b Galactic photoelectric 
absorption extracted from the HEASARCH Database (LAB survey Kalberla et al. 2005); 



c net exposure times for the three EPIC instruments after data cleaning in ks; 
d extraction radius is half the cooling radius from Peres et al. (1998); c the two 
observations of Abell 4059 were merged together before the spectral analysis. 



is the radius encompassing a spherical density contrast of 500 
with respect to the critical density. By converting our extraction 
radii in 1-500 units we find a rather peaked distribution centered 
around 0.08 rsoo- 



4. Spectral analysis 

In this paper we wish to measure abundances for Si, Fe and Ni. 
To achieve robust estimates we choose to make our measures 
from the K-shell lines only, avoiding L-shell blends where uncer- 
tainties associated to the atomic physics are larger. We achieve 
this by setting a lower energy threshold of 1.8 keV which in- 
cludes Si K-shell emission (2 keV) and excludes the Fe and Ni 
L-shell blends (~ 1.-1.5 keV). 

The background subtraction for the high surface brightness 
core regions of cool-core clusters is less critical than in the case 
of more external low surface brightness regions. We therefore 
have decided to subtract the background using blank-sky fields 
instead of proceeding with a more detailed modeling of the dif- 
ferent background components. The blank-sky fields f or EPIC 
MOS and pn were produced by iLeccardi & Moleridil (|2008a) 



(see Appendix B "The analysis of blank field observations" in 
their paper) by analysing a large number of observations for a 
total exposure time of ~ 700 ks for MOS and ~ 500 ks for pn. 

We have refined our background analysis performing also a 
background rescaling for each observation separately to account 
for temporal variations of the background. We have estimated 
the background intensity from a spectra extracted from an ex- 
ternal ring between 10' and 12' centered on the emission peak, 
taking into account only the 10-12 keV band (to avoid possi- 
ble extended cluster emission residuals in this region). We have 
than rescaled the blank-sky fields background to the local value. 
This procedure is important especially when deriving the nickel 
abundance from its emission lines at ~ 8 keV. Indeed, in this 
hard spectral region both the EPIC effective area and the sur- 
face brightness of the relatively low temperature cluster cores de- 
crease rapidly, and, the background becomes progressively more 
important. 

All the spec tral fits were p erformed with the XSPEC package 
(version 11. 3.2. lArnaudlll996l) . 

We have analyzed each cluster spectra with three differ- 
ent models: (1) a one temperature thermal model with the 
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6.3 6.4 6.5 6.6 6.7 6.8 6.3 6.4 6.5 6.6 6.7 6.8 




channel energy (keV) channel energy (keV) 



Fig. 1. EPIC-pn spectral data (crosses) and best-fit models (continuum lines) for 3 high counting statistics clusters: Perseus (on the 
top), 2A 0335 + 096 (in the middle) and Centaurus (at the bottom). Panels on the left show results from a 4T model fit without 
Gaussian smoothing whereas plots on the right are smoothed. In all cases we also plot the data divided by the folded model. The 
spectra are rebinned to allow graphical clarity. The energy range of each plot is 600 eV centered on the 6.7 keV line (redshifted at 
the cluster distance). In the panels on the left are clearly visible the "wings" around the Fe line in the residuals that are an indication 
of the resolution miss-calibration of the redistribution matrix of the pn (see discussion in Sect. 4.1). 



plasma in collisional ionization equilibrium (vtnekal model 
in the XSPEC nomenclature), referred as IT model there- 
after, (2) a two temperature thermal model (vmekal+vmekal), 
2T model thereafter, and, (3) a multi-temperature model 
(ymekal+vmekal+vmekal+vmekal), 4T model thereafter. 

All models have been multiplied by the Gal actic hydrogen 
colum n density, Ng, determined by HI surveys dKalberla et al.l 
120051) through the wabs absorption model in XSPEC. In the fol- 
lowing sections we report results from spectral fits with Nh fixed 
at the weighted Galactic value only. We have also allowed Nh to 
vary in all three models finding no significant differences in the 
derived silicon, iron and nickel abundance values. 

The redshifts have been taken from the NASA/IPAC 
Extragalactic Database and have been always left as free param- 
eters to account for small calibration differences between the two 
MOS and the pn. We have allowed the temperatures of IT and 
2T models to vary freely, whereas in the 4T model temperatures 
of the four components were fixed at 1, 2, 4 and 8 keV, respec- 
tively. The normalizations of the models and the elements Si, S, 
Fe and Ni were left free to vary. Ar and C were jointed to S as 



these elements are less abundant and we do not plan to study 
them in the course of this work. 

As pointed out bv lLeccardi & Molendil(l2008al) . when fitting 
the spectra with XSPEC it is opportune to allow the metallicities 
to assume negative values. This procedure is necessary to avoid 
underestimates on the derived metal abundances which could af- 
fect the measurements especially for the case of low metallicity, 
statistically poor spectra (for a more detailed discus sion of this 
point see Appendix A in lLeccardi & Molen di 2008a). 

Abund ances are measured relati ve to the solar photospheric 
values of lAnders & Grevessd (1 19891) . where Fe = 4.68 x 10~ 5 , 
Si = 3.55 x 10~ 5 and Ni = 1.78 x 10~ 6 (by number relative to 
H). We have chosen these values to allow direct comparison with 
other works present in the literature. Spectra from all three EPIC 
instruments (MOS1, MOS2 and pn) were fit individually. 

In the course of our analysis we have found that spectra with 
less than about 3500 source counts could not be used to constrain 
the abundances of nickel and silicon, we therefore have decided 
to eliminate from the sample all clusters having less than 3500 
source counts in the adopted energy band (i.e. 1.8-10 keV) 
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Table 2. Relative differences between metal measurements made with the three EPIC instruments MOS1, MOS2 and pn computed 
over the whole sample. For each element in the first line only statistical errors are accounted for, in the second line a 3% systematic 
error is summed up in quadrature to each abundance measurement (see discussion in Sect. 4.3). 



mekal model 


Metal X 


(XmOSI ~ X M os2)/X M os2 


(Xmosi ~ X p „)/X p „ 


(X-MOS2 ' 


X pn) 1 pn 


Fe 


0.02 


+ 0.01 


0.07 ± 0.01 


0.02 


±0.02 


+3% syst.err. 


0.02 


±0.02 


0.04 ± 0.02 


0.01 


±0.02 


Si 


0.01 


+ 0.02 


0.19 ±0.03 


0.10 


±0.03 


+3% syst.err. 


-0.03 


±0.04 


0.11 ±0.04 


0.08 


±0.04 


Ni 


-0.12 


±0.08 


0.14 + 0.11 


0.01 


±0.09 


+3% syst.err. 


-0.13 


±0.08 


0.12 ±0.11 


0.00 


±0.10 



apec model 


Metal X 


(Xmosi - Xmos2)/Xmos2 


(X-MOSl ' 


" X p „)/X p „ 


(X-MOS2 - 


*^pn)/*^pn 


Fe 


0.03 


±0.01 


0.07 


±0.01 


0.02 


±0.02 


+3% syst.err. 


0.02 


±0.02 


0.04 


±0.02 


0.01 


±0.02 


Si 


0.01 


±0.03 


0.22 


±0.03 


0.13 


±0.04 


+3% syst.err. 


-0.03 


±0.05 


0.13 


±0.04 


0.11 


±0.04 


Ni 


-0.16 


±0.09 


0.13 


±0.14 


0.01 


±0.11 


+3% syst.err. 


-0.17 


±0.10 


0.12 


±0.14 


0.01 


±0.11 



in at least one of the three EPIC instruments. The five clusters 
removed from the sample by adopting this criterion are: Abell 
1644, Abell 1651, Abell 2065, Abell 2063 and Abell 3562. 

We have also excluded Cygnus A because of its peculiar 
core emission. A powerful double-lobed radio galaxies, QSO 
B1957+405 featuring huge jets feeding very large hotspot re- 
gions which are well detected in X-rays, resides in the core of 
this cluster. In the light of the difficulty of removing efficiently 
the emission associated to the radio galaxy from the spectra of 
the core region we have preferred to exclude this cluster from 
the sample. 

In the case of the nickel abundance measurement we have 
introduced a further selection criterion to exclude clusters which 
are background dominated in the hard part of the spectrum. We 
have selected the 7-9 keV energy band, redshifted at the clusters 
distance, around the nickel line at ~ 8 keV, and we have mea- 
sured the relative difference between the source and background 
count rates, (cr sou - crh g k)lcrbk g , in that band. We have then con- 
sidered only nickel measurements from clusters with relative dif- 
ference larger than 1, namely clusters with source count rates 
which were at least twice the background count rates in that hard 
band. Applying this selections Virgo, Abell 262 and Abell 1060 
are also excluded from the sample. It is worth noting that our 
pn spectra are only mildly contami nated by the fluore scence line 
complex located around 8 keV (e.g. lNess et alj2009l) . the reason 
being that our spectra are extracted from the central regions of 
the pn detector where, due to a hole in the pn electronics box, 
the fluorescence lines are weak. 

In summary we have measured Si and Fe abundances from a 
sample of 26 clusters and Ni from a subsample of 23 systems. 

4. 1 . Addendum to the spectral analysis of the pn 

iMolendi & Gastaldello ( 2009) analyzing the pn spectra of a long 
XMM-Newton observation of the Perseus cluster first noted in the 
best-fit residuals the presence of a substantial structure around 
the Fe Kor line (see Sect. 3.1.1 and Fig. 2 in their paper). This 
structure is attributed to an incorrect modeling of the pn spectral 
resolution within the redistribution matrix. 



Investigating this features extensively in the pn spectra of our 
cluster sample we have found that the resolution miss-calibration 
can be compensated for by including a multiplicative compo- 
nent that performs a gaussian smoothing of the spectral model 
(gsmooth in XSPEC). We set the width of the gaussian kernel to 
be 4 eV (FWHM) at 6 keV and assume a power-law dependency 
of the width on the energy with an index of 1 . 

Figure 1 shows the pn spectra and best-fit residuals for three 
bright clusters: Perseus, 2A 0335+096 and Centaurus clusters. 
We have fitted the spectra with the 4T model without and with 
the gaussian smoothing. In the former cases (left panels in the 
Figure) prominent residuals are present in the "wings" of the Fe 
line, while in the latter (right panels) the residuals are signifi- 
cantly reduced. In 13 (11) out of 26 cases the modified mod- 
els applied to pn spectra provide a substantial better fit (i.e. 
Ax 2 > 2.7(4.0)) than the un-modified models. We note that Si 
and Ni lines do not show similar residuals with respect to the 
best fit model, most likely because the statistical quality of the 
data is not as high (Si and Ni) and because the spectral resolution 
is significantly poorer (Si). 

After the application of the gaussian smoothing the cross- 
calibration between the pn and MOS improve. For instance, 
the relative differences, computed over the whole sample, be- 
tween iron abundances estimated from pn and MOS1 spectra, 
i.e. (FeMosi - Fe / ,„)/Fe p „, with the un-smoothed and smoothed 
4T model decrease from 12% to 7%, whereas the ones between 
pn and MOS2 decrease from 9% to 4% (errors on the given per- 
centages are always 1%). 

4.2. Choice of the best spectral model for each cluster 

In this section we discuss how we have selected the model that 
best fits the spectral data among the IT, 2T and 4T models de- 
scribed above. 

Applying the statistical F-test between IT and 2T (or IT and 
4T) models, we find that only the first 6 out of 26 clusters with 
larger statistics (i.e. Perseus, Virgo, 2A 0335 + 96, Centaurus, 
Abell 478 and Abell 4038) show overwhelming evidence of 
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multi-temperature structure, namely an F-test probability, PF, 
smaller than 1% for all EPIC instruments. 

The majority of the other clusters (17/26 for 2T model and 
22/26 for 4T model) however, show a significant improvement 
(PF < 1%) in at least one of the EPIC instruments. For these 
clusters the relative differences between the metal abundances 
measured in each instrument with the IT and 2T (or IT and 4T) 
models are always smaller than 2-3 cr for Fe and Si, and, consis- 
tent with zero for Ni. Therefore, from a purely statistical point of 
view, the choice of IT or multi-temperature (2T or 4T) models 
results in modest differences. Nevertheless, we point out that our 
spectra are extracted from the very central regions of cool-core 
clusters a nd that these regions always display temperature gradi- 
ents (e.g. iLeccardi & Molend i 2008b and references therein). It 
follows that the best description of the spectrum of this plasma 
is through a multi-temperature model. For this reason and for 
the F-test results above, we decide to use a multi-temperature 
modeling for all the clusters in the sample. 

An F-test between 2T and 4T models is not possible as 
the two models have the same numbers of degrees of freedom. 
However, relative differences between the metal abundances 
measured with the two models show that the systematic uncer- 
tainties associated to the different modeling of the data are below 
2-3% for Fe and Si and consistent with zero for Ni. These sys- 
tematic errors are of the same order or smaller than systematic 
errors on the abundances due to calibration differences between 
the three EPIC detectors (systematic errors given by the cross 
correlation between EPIC instruments will be discussed here- 
inafter in Sect. 4.3). Therefore the two models are substantially 
equivalent to describe our data. We choose to use for all the clus- 
ters the 4T model. 



Table 3. Measurements o f silicon, iron and nickel a bundances 
relative to the solar value (lAnders & Grevesse1 IT989) for the fi- 
nal sample (26 clusters for Si and Fe, 23 for Ni) from the 4T 
mekal spectral model. Clusters are sorted by decrescent counts 
statistics in the pn. Errors on abundances include a 3% system- 
atic error as detailed in Sect. 4.3. 
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12 


Ahell 71 QQ 


0.92 


±0.10 


0.45 


±0.02 


1.11 


±0.41 


1 T. 
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Ahell 707Q 


0.61 


±0.15 


0.58 


±0.03 


1.66 


±0.86 


14 


Abell3112 


0.77 


±0.10 


0.58 


±0.04 


1.17 


±0.43 


15 


HYDRA A 


0.46 


±0.11 


0.39 


±0.02 


0.75 


±0.13 


16 


Abell 496 


0.85 


±0.06 


0.55 


±0.02 


1.72 


±0.37 


17 


AWM7 


1.13 


±0.07 


0.64 


±0.02 


1.24 


±0.30 


18 


ABELL 4059 


0.78 


±0.03 


0.59 


±0.04 


1.00 


±0.12 


19 


Abell 3571 


0.52 


±0.15 


0.48 


±0.03 


1.81 


±0.09 


20 


Abell 1650 


0.36 


±0.16 


0.57 


±0.03 


1.11 


±0.62 


21 


Abell 1689 


0.40 


± 1.68 


0.33 


±0.03 


1.96 


±0.33 


22 


PKS 0745 - 191 


0.61 


±0.22 


0.43 


±0.03 


0.49 


±0.22 


23 


Abell 262 


0.93 


±0.06 


0.53 


±0.06 






24 


Abell 2204 


0.83 


±0.24 


0.53 


±0.03 


1.27 


±0.49 


25 


Abell 85 


0.94 


±0.23 


0.55 


±0.02 


2.49 


±0.74 


26 


Abell 3558 


0.84 


±0.36 


0.49 ± 0.02 


1.68 


± 1.15 



4.3. EPIC Cross-calibration Issues 

For each cluster we compute the relative differences between 
iron, silicon, and nickel measured with MOS1, MOS2 and pn, 
and subsequently compute the weighted averages of these dif- 
ferences for the whole sample (i.e. 26 clusters for Fe and Si, 
and 23 clusters for Ni, results are from 4T model). The result- 
ing mean values for Si, Fe and Ni are shown in Table [2] From 
the results shown in this table it is clear that systematic errors 
are often dominant with respect to statistical ones for Fe and Si, 
whereas Ni appears to be fully dominated by statistical uncer- 
tainties. Another important point is that the relative differences 
between MOS 1 and MOS2 are smaller with respect to the differ- 
ences between MOS1 or MOS2 and pn. 

To account for systematic differences between measures ob- 
tained from different detectors and spectral models we sum in 
quadrature a 3% systematic error to each Fe and Si abundance 
measurement. The relative differences between iron, silicon, and 
nickel measured with MOS1, MOS2 and pn computed over the 
whole sample, including the 3% systematic errors, are smaller 
than those computed without them and always significant at less 
than 3<r (these values are given in the second line relative to each 
element in Table |2j. 

For each cluster we compute an EPIC Fe and Si abundance 
by performing error weighted averages over the three detectors, 
with errors including the 3% systematic described above. Errors 
on the EPIC Fe and Si abundances are computed by dividing 
the error weighted standard deviation of the EPIC abundance by 
the square root of 3. To avoid errors on the EPIC abundances 
from falling below the systematic level a 3% systematic error is 
summed in quadrature. 



As shown in Table [2] nickel measurements are clearly dom- 
inated by statistical errors, nevertheless we decided for internal 
consistency to compute the nickel value for each cluster as de- 
scribed above for Si and Fe. Indeed we expect Ni to be prone to 
the same systematics affecting Si and Fe. 

The distribution of Fe and Si abundances is not very sym- 
metric, this entails that it is not particularly well represented 
by a gaussian. We have determined that the major cause is the 
presence of a few measurements with extremely small errors. 
Introducing a systematic error of 3% (see above) alleviates the 
problem considerably. The fact that the data is not normally dis- 
tributed implies that quantities such as the mean and standard de- 
viation do not enjoy the properties they "normally" do. A man- 
ifestation of non-normal behavior can be observed in Table 2: 
consider for example the Fe abundance measured with mekal, al- 
though the mean relative difference for MOS 1 and MOS2 mea- 
sures is small 2%, and the mean relative difference for MOS1 
and pn is large, e.g. 7%, the mean relative difference for MOS2 
and pn is modest, e.g. 2%. Similar results are observed for Si 
measured with mekal and with Fe and Si measured with apec. 

The final abundances are reported in Table [3] column (1) 
is the cluster name, column(2), (3) and (4) are the EPIC error 
weighted averages of Si, Fe and Ni, respectively, with their l<x 
errors. 

For completeness we have also computed all the abundances 
presented in Table [3] using the 2T model, as expected we found 
only modest differences with respect to those estimated with the 
4T model. The relative differences between the mean Fe, Si and 
Ni values computed for the whole sample with the 2T and 4T 
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models are -3% + 1%, 2% + 7% and -1% + 9%, respectively. 
In the case of Fe we detect a systematic which is of the same 
order of the one associated to the choice of detector, while in the 
case of Si and Ni indeterminations are sufficiently large to hide 
systematics of the same order. 

4.4. apec vs. mekal 

Mekal is not the only collisional ionization equilibrium (CIE) 
emission model in XSPEC. For quite some tim e an alternative 
plasma code known as apec (iSmifh et al.ll200ll) has been avail- 
able. Some authors have performed comparis ons between spec- 
tral fits performed with the two codes (e.g. ISanders & Fabian 
2006; Ide Plaa et al.ll2007l) . I Sanders & Fabianl (120061) . analyzing 
a long Chandra observation of the core of Centaurus find that 
while there are little differences between Si and Fe abundances 
estimated with the two codes, the Ni abundance found with apec 
is somewhat smaller than that found with mekal, a similar result, 
albeit with much less statistics, is found by de Plaa and collabo- 
rators. Here we wish to investigate how the differences between 
the plasma codes impact on our abundance and abundance ra- 
tio measurements. To this end, we reanalyze all our objects with 
a multi-temperature model (vapec+vapec+vapec+vapec) analo- 
gous to the 4T mekal model which we dub 4T apec. A compar- 
ison amongst apec metal abundances measured with the three 
EPIC instruments (see Tab.[2]i provides results that are very sim- 
ilar to those found for mekal; including a systematic of 3% on 
the measures alleviates difference between the three instruments 
in much the same way (see Tab.|2]i it does for the 4T mekal fits. 
Mean metal abundances for Si, Fe and Ni and abundance ra- 
tios, Si/Fe and Ni/Fe obtained with apec and associated intrinsic 
scatters are reported in Tab. |4] Comparing the sample averaged 
metal abundances measured with mekal with those obtained with 
apec we find that: Fe is almost unchanged, Fe a pecfPe m ekai = 
1.05 + 0.01; Si is somewhat higher, Si apec /Si me k a i = 1.11 ± 0.02 
and Ni is lower, ^i ape c^mekai — 0.82 + 0.04. The Si/Fe ratio, as 
measured with apec, is slightly larger than that estimated with 
mekal, (Si/Fe) flpe , c /(Si/Fe) me , J ( :a / = 1 .06 + 0.02, while Ni/Fe is sub- 
stantially smaller, (Ni/Fe) apec /(Ni/Fe) mefa i = 0.77 + 0.04. 

In summary, while for Fe and Si/Fe the systematics associ- 
ated to the cross-calibration between EPIC experiments, 3%, are 
comparable to those related to the emission model, 5% and 6% 
respectively, for Si, Ni, and Ni/Fe the dominant source of error 
are the systematics related to the emission model, 11%, 18% and 
23% respectively. It goes almost without saying that the statis- 
tical errors on the mean Si, Fe, Ni, Si/Fe and Ni/Fe are always 
dominated by systematics. 

5. Abundances results in Cool-Cores 

5. 1 . Mean Iron, Silicon and Nickel Abundances 

In Fig. |2]we show the abundances of silicon (upper panel), iron 
(middle panel) and nickel (lower panel) in the cool-cores of our 
clusters (we use hereinafter results from the mekal model al- 
though results from both spectral codes are reported in Tab. 3J. 
Visual inspection of the plots indicates that the abundance distri- 
butions show only small, 20%-30% variations with respect to the 
mean. The only object that is clearly out of the distributions is 
the Centaurus cluster whose core is uncommonly rich in metals 
(e.g. lSanders & Fabianll2006l) . 

To quantify the qualitative indications shown in Fig. [2] we 
have estimated the properties of the abundance distributions with 
a maximum likelihood algorithm that postulates a gaussian par- 



Table 4. Averages of Si, Fe, Ni abundances and Si/Fe, Ni/Fe ra- 
tios for the cool-core cluster sample in solar units. Mean values 
are given both for the total sample and with the exclusion of the 
peculiar Centaurus cluster. Average abundance values and intrin- 
sic scatters are computed with the maximum likelihood method 
described in lMaccacaro et alj dl988l) . 



mekal model 


Metal 


mean 


scatter 


Si 


0.75 


±0.03 


0.22 


±0.03 


Si no Cent. 


0.72 


±0.03 


0.17 


±0.02 


Fe 


0.51 


±0.02 


0.11 


±0.01 


Fe no Cent. 


0.49 


±0.01 


0.07 


±0.01 


Ni 


1.28 


±0.08 


0.45 


±0.06 


Ni no Cent. 


1.18 


±0.07 


0.34 


±0.06 


Si/Fe 


1.47 


±0.05 


0.27 


±0.05 


Si/Fe no Cent. 


1.45 


±0.05 


0.30 


±0.05 


Ni/Fe 


2.41 


±0.13 


0.60 


±0.12 


Ni/Fe no Cent. 


2.40 


±0.14 


0.63 


±0.13 




apec model 


Metal 


mean 


scatter 


Si 


0.83 


±0.04 


0.23 


±0.03 


Si no Cent. 


0.80 


±0.03 


0.19 


±0.03 


Fe 


0.53 


±0.02 


0.12 


±0.01 


Fe no Cent. 


0.51 


±0.01 


0.08 


±0.01 


Ni 


1.07 


±0.10 


0.53 


±0.07 


Ni no Cent. 


0.96 


±0.08 


0.39 


±0.06 


Si/Fe 


1.56 


±0.05 


0.25 


±0.05 


Si/Fe no Cent. 


1.56 


±0.05 


0.26 


±0.05 


Ni/Fe 


1.88 


±0.14 


0.67 


±0.14 


Ni/Fe no Cent. 


1.83 


±0.15 


0.67 


±0.15 



ent distribution describ ed by a mean and an intrinsic dispersion 
(Maccaca ro et alll98 8). Our analysis shows that our abundance 
distributions feature a substantial intrinsic dispersion, roughly 
30% for Si and Ni and 20% for Fe. Best fitting means and in- 
trinsic dispersions for all species are summarized in Table [4] 
To examine the influence of the outlier Centaurus cluster on the 
averages and scatters, we have redone the computation exclud- 
ing it from the sample, results are reported in Tab. [4] We do not 
find any dramatic change in the derived quantities, especially in 
the scatters of the abundance distribution that are only slightly 
smaller than those computed with the Centaurus cluster. 

We have compared ou r Fe abundance estimates with those 
from lLeccardi et al] (120091) who made measures in a circular re- 
gion with radius 0.05 rise, which roughly corresponds to our 1/2 
icooi extraction radius. We have considered their cool-cores sys- 
tems excluding objects contained in our sample (these selections 
leads to a subsample of 13 clusters located in the 0.1-0.25 red- 
shift interval). Applying the maximum likelihood algorithm to 
their Fe distribution we find a mean Fe abundance of 0.49 + 0.02 
and a relative intrinsic scatter of 16% + 3%, which are both in 
excellent agreement with our measurements. 

The mean Si and Fe abundance values r eported in T ableware 
also in good agreement with what found by Tam ura et alj {2004) 
for the innermost regions of their cool-c ore clusters with tem per- 
atures larger than 3 keV (see Table 6 inlTamura et al.l l2004). 

On the contrary Ide Plaa et all d2007l) found Si, Fe and Ni 
abundances somewhat smaller than ours. By applying the max- 
imum li kelihood method to the abundances given in their Table 
A. 1 (see lde Plaa et al.l2007l) . we find an average Si of 0.39+0.02, 
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Fig. 2. Silicon (upper panel), iron (medium panel) and nickel (lower panel) abundances measured in the core region with r < 0.5 
r C ooi for the cluster sample. The total number of clusters is 26 for Si and Fe, and 23 for Ni. Clusters are sorted by decreasing pn 
counts statistics according to Table [3] Weighted mean values and intrinsic scatters of the relations are given in Table |4] 



Fe of 0.35 ± 0.02 and Ni of 0.76 ± 0.05 (here and in the rest of 
this work abundances tak en from the literature are rescaled to 
lAnders & Grev esse 1989), with intrinsic scatters for Si and Fe 
of 24% and 17%, respectively. The smaller abunda nce values are 
not su rprising since the extraction regions used by de Pla a et all 
J2007h are about three times larger than those considered in our 
work and because their sample contains both cool-core and non 
cool-core clusters. 

Interestingly, the mean Si and Fe found by 
iRasmussen & Ponmanl (120071) from a sample of galaxy 
groups are very similar to ours. From their Fig. 1 1 (panels b and 
c) we estimate a mean Si of ~ 0.75 and a mean Fe of ~ 0.50 
within 0.08 r/rsoo which roughly corresponds to our 1/2 r coo i 
extraction radius. 



5.2. Silicon-to-lron and Nickel-to-lron ratios 

In Fig. [3] we plot the silicon-to-iron (upper panel) and nickel-to- 
iron (lower panel) ratios distributions for our objects. Using the 
maximum likelihood algorithm described in the previous section 
we have computed the mean and the intrinsic dispersion for the 
Si/Fe and Ni/Fe distributions, results are shown in Tab. [4] The 
relative scatter found for Si/Fe is about 20%, i.e. comparable to 
the one estimated for Fe and somewhat smaller than the one for 
Si; for Ni/Fe the scatter is only slightly larger 25%. Interestingly 
the Centaurus cluster, which features very high central abun- 
dances, is characterized by Si/Fe and Ni/Fe ratios very similar 
to those observed in other clusters. 

As expected from the results on the Si and Fe abundances 
shown in the previous Section we find that our ave rage Si/Fe ra- 
tio is in good agreement with the value found by Tamur aet al.l 
(2004), 0.16 + 0.2, in the innermost bin of their high temperature 



(> 3 k eV) cool-core clusters sample, whereas in Ide Plaa et al.l 
^20071) the mean Si/Fe and Ni/Fe abundance ratios are 1 .07+0.03 
and 2.07 + 0.13, which are significantly lower than our values 
(the abundance ratios are computed with the ma ximum likeli- 
hood method applied to the data given in Tab A. 1 o flde Plaa et al.l 
120071) . As already mentioned, the obs erved difference be tween 
our abundance ratios and those from de Plaa et a il (120071) could 
derive from the different sample selection. 



We have compared our Si/Fe ratio with those m easured 
for a sample of groups (Rasmus sen & Pon man 200 7)) and o f 
X-ray luminous elliptical galaxies ( Humphrey & Buotel 12006). 
We note that in both papers abun dances are reported in so- 
lar units which differ from ours, iGrevesse & Sauvall J1998J) 
for |R asmussen & Ponman| (l2007l) and lAsplund et alj d2005l) for 
Humphrev & Buote (2006); to allow immediate comparison 
w ith our results the values s hown below have been converted 
to lAnders & Grevessd ([1989) units. The measure on the groups 
at scales comparable to thos e adopted for our clusters, 0. 1 rsoo, 
(see Fig. 1 1 f and Sect. 5.3 in Rasmussen & Ponmanll2007l) . pro- 
vides a mean Si/Fe ratio of 1.35 and a relative scatter of 32% 
which are both in broad agreement with ours. The measure on 
the galaxies was made us ing data from high-luminos ity ellipti- 
cal reported in Table 2 by Humphrey & Buotel (|2006). We have 
applied the maximum likelihood method used for our own data 
deriving a mean Si/Fe value of 1.50 + 0.07 and a relative intrinsic 
scatter of 16% ± 3%, which are both in good agreement with our 
estimates for cool-cores. 



Summarizing, the average Si/Fe ratio appears to be nearly 
constant from the galactic to the rich clusters scale suggesting a 
common enrichment scenario in these objects. 
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Fig. 3. Silicon-to-iron (upper panel) and nickel-to-iron (lower panel) ratios for the cluster sample. 



6. SNIa vs. SNcc 

There are indications that both type la and core-collapsed 
supernovae contribute significantly to the enrichment of 
the i ntergalactic medium in co re s of cool-core c lusters 
(e.g. Ilshimaru & Arimotol 119971 : iFinoguenov et alj [2000; 
iGastaldello & Molendi |2002|) . We try here to determine the 
relative proportion between the two SNe types using the 
silicon-to-iron and nickel-to iron abundance ratios normalized 
to the solar value. 

The observed X,/Fe ratio, where X, is the i-th element, 
can be expressed as a linear combination of the X,/Fe ra- 
tio expected from type la an d core-collapsed supernovae (e.g. 
Gastaldello & Molendi 2002): 



XA 

Fe /obs 



(!) 



♦Ci-fHg) • 

SNIa \ lc /SNcc 



where 



M, 



Fe.SNIa 



M, 



Fe.SNIa 



+ M, 



Fe.SNcc 



(1) 



(2) 



is the SNIa iron-mass-fraction, i.e. the fraction of the iron 
mass synthesized in type la supernovae, M Fe ,sNia an d -^Fe.SNcc 
are respe ctively the total Fe mass ejected by SNIa and SNcc. 

iMatteucci & Chiappinil 12005) caution against an inappropri- 
ate use of equation (1) as in principle this equation is valid only 
if one is interested in the global Fe production (i.e., Fe in stars, 
galaxies and the ICM). Indeed the model leading to eq. (1) does 
not consider the finite timescale over which the gas is processed 
through stars and ignores the mechanism of chemical enrich- 
ment of the ICM fr om galaxies. However, a s alrea dy noted by 
Ide Plaa et ail d2007l) and lHumphrev & Buotd d2006l) . eq. (1) can 
be applied to the ICM, quite simply the iron mass fractions that 
are derived from it should be interpreted as the iron mass frac- 



tions that would be needed to enrich the ICM, not as the actual 
iron mass fraction produced throughout the history of the cluster. 

By inserting Si/Fe and Ni/Fe ratios expected for SNIa and 
SNcc and the mean observed Si/Fe and Ni/Fe ratios in eq. 
(1) we estimate the SNIa iron-mass-fraction We compute 
£ by combining several theoretical SNe yields, the ones for 
SNIa origi nate from two physically different sets of models 
taken from Iwamot oet alJ (1999), namely the slow deflagration 
(W7 and W70) and the delayed detonation (WDDs and CDDs) 
explosion models. For core-col l apsed supernovae we us e the 
yields fromlNomoto et aD(l2006h . lChiem & Limongil d2004 and 
IWoosley & Weaverl ri 19951) . details on the models and the com- 
putation of the yields are provided in App. A. Since we have 
two equations, one for Si/Fe the other for Ni/Fe, with one un- 
known, we are in a position to determine whether a given 
combination of la and core-collapsed models can adequately re- 
produce the observ ed ratios. One way of going about this (e.g. 
Ide Plaa et alJ2007T) is to perform a^- 2 fit using the observed Si/Fe 
and Ni/Fe ratios and associated errors as the data, the Si and Ni 
yields relative to Fe for SNIa and SNcc, namely (Si/Fe)sNia and 
(Si/Fe)sNcc, as constants and £ as fitting parameter. We prefer 
to proceed differently, indeed it has been noted that the uncer- 
tainties associat ed to the SNe yields are i n the order of tens 
of percent (e.g.. IWoosley & Weaverl 1 1995b iGibson etalJll997l 
P. Young priv. comm., and the scatter in the yields we have 
gathered from the literature fully confirms this, see Table I A. 1 1 
and Appendix A), while our observed mean values of Si/Fe 
and Ni/Fe have errors smaller than 5%, if we, for the time be- 
ing, neglect the uncertainties associated to the spectral model. 
Consequently a robust estimate of £ should first of all take into 
account the uncertainties in the yields. Given the lack of pre- 
cise information we have assumed a uniform indetermination in 
the Si and Ni yields relative to Fe for SNIa and SNcc of 20%. 
For each combination of SNIa and SNcc models we allow the Si 
and Ni yields relative to Fe for SNIa and SNcc to vary by 20% 
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Table 5. Supernovae type la mass fraction, £ = Mp e £mal{MFe,SNia + Mf ei sNcc), obtained for different theoretical SNe silicon-to- 
iron and nickel-to-iron ratios from the observed Si/Fe and Ni/Fe ratios. SNI a yields are computed from the full set of models for 
deflagration (W7, W70) and delayed detonation (WDDs, CDDs) scenarios by Iwamoto et al.l d 19991) . SNcc yields are integrated over 
a Salpeter (x = 1.35) or a A Y top-heavy (x = 0-95) IM F, with abun dance of the progenitor Z =0.02. Models for SNcc are taken from 
iNomoto et al. 1 (120061 No061. IChieffi & Limongil (2004. CL04) and lWooslev & Weaverl(ll995l W951. 



mekal model 


SNIa 


SNcc 


W7 


W70 


WDD1 


WDD2 


WDD3 


CDD1 


CDD2 


No06,S 


0.528 -0.571 


0.515-0.734 












No06,T 


0.575 - 0.579 


0.562 - 0.760 












CL04,S 




0.505 - 0.582 




0.546 - 0.789 


0.505 - 0.728 




0.534 - 0.770 


CL04,T 




0.559 - 0.728 




0.600 - 0.727 


0.559 - 0.758 




0.588 - 0.727 


W95,S 






0.616-0.677 


0.518-0.775 


0.476-0.712 


0.635 - 0.652 


0.505 - 0.756 


W95,T 








0.554-0.792 


0.512-0.732 




0.541 -0.774 



apec model 


SNIa 


SNcc 


W7 W70 


WDD1 


WDD2 


WDD3 


CDD1 


CDD2 


No06,S 


0.475 -0.511 




0.514-0.764 


0.475 - 0.707 




0.502 - 0.747 


No06,T 






0.564 - 0.789 


0.526 - 0.735 




0.553 -0.773 


CL04,S 




0.591 -0.886 


0.502 - 0.759 


0.581 -0.701 


0.608 - 0.841 


0.490 - 0.741 


CL04,T 




0.647 - 0.855 


0.561 -0.787 


0.523 - 0.734 


0.663 - 0.799 


0.549-0.771 


W95,S 




0.560-0.885 


0.640 - 0.744 




0.577 - 0.885 


0.640 - 0.725 


W95,T 




0.598 -0.895 


0.561 -0.763 




0.615-0.868 


0.561 -0.745 



and, using the observed Si/Fe and Ni/Fe ratios, we determine the 
range of values of £ for which eq. (1) is satisfied. For some com- 
binations of SNIa and SNcc models, despite the generous 20% 
range, there are no values of £ satisfying eq. (1), which implies 
that the observed Si/Fe and Ni/Fe cannot be reproduced for that 
particular combination of SNIa and SNcc models. 

In Table [5] we report ranges of values for £ for all possible 
combinations of la and core-collapsed supernovae models. The 
overall permitted range for £ is large 0.48 - 0.79 and sensitive 
to the indetermination in the SNe yield ratios, indeed assuming 
10% or 30% rather than 20% we derive a range of 0.55 - 0.73 
and 0.37 - 0.85, respectively. Since the indetermination in the 
abundance ratios related to the emission model is comparable to 
the one in the yields we recompute the SNIa iron-mass-fraction 
assuming abundance ratios derived from apec rather than mekal. 
Results are rather similar, see Table IA.11 the overall permitted 
range for £ is similar 0.49 - 0.90, the major difference is in which 
of the combination of SNIa and SNcc models provide valid so- 
lutions. 

Equations (1) and (2) may be generalized to any atomic 
species and can therefore be used to obtain the X, elements SNIa 
gas-mass-fraction. We have done so and derived the Si and Ni 
SNIa gas-mass-fraction ranges which are respectively 0. 14-0.49 
and 0.26 - 0.88. 

By adopting abundance ratios derived from apec we find 
broader ranges both for Si, 0.12- 0.86, and Ni, 0.01 - 0.90. The 
reason for this difference rests in the smaller Ni/Fe ratio found 
with apec. More specifically the mekal value, 2.41, is relatively 
high when compared to the Ni/Fe ratios predicted by SNIa and 
SNcc models. This implies that the observed Ni/Fe and Si/Fe ra- 
tios may only be reproduced by a limited combination of SNIa 
and SNcc models. Conversely the apec value, 1.88, is closer to 
the mean Ni/Fe ratio predicted by models and may be repro- 
duced by a broader combination of SNIa and SNcc yields. 

The total Fe mass ejected by SNIa, Mpe.sNia, and by SNcc, 
A^Fe.SNcc can be rewritten as: 

A^Fe.SNIa = A^SNIa jFe.SNIa (3) 



and 

A^Fe.SNcc = A^SNcc jFe.SNcc, (4) 

where A^sNia and A^sncc are respectively the total number of 
type la and core-collapsed supernovae and yFe.SNia and jFe.SNcc 
are respectively the Fe yields per SNIa and SNcc. 

We define, / = -/VsNia/O^VsNia + A^sncc) as the number ratio of 
supernovae type la over the total number of supernovae; com- 
bining eq. (2) with eqs. (3) and (4) we solve for /: 

j, _ £, jFe.SNcc ^ 

£, jFe.SNcc + (1 - jFe.SNIa 

In Table [6] we report estimates for / obtained by inserting 
allowed ranges for £ in eq. (5). We do not introduce any fur- 
ther indetermination in the Fe yields as the scatter in Si and Ni 
yields is much larger than the one in Fe and already accounted 
for through the scatter introduced in the Si and Ni yields relative 
to Fe. As for £ the estimates are expressed in the form of per- 
mitted ranges. The overall permitted range for / is rather large: 
0.10-0.38. Using abundance ratios estimated with the apec code 
we estimate a somewhat broader permitted range of 0.09-0.61. 

7. Discussion 

In this work we have analyzed the central regions (r < r coo i/2) of 
a sample of 26 nearby and rich cool-core cluster s derived from 
the original B55 cluster sample (Edge et al. 1990). We have used 
these regions characterized by the highest cluster photon statis- 
tics to identify all the possible sources of systematic uncertain- 
ties in the abundance measurements of the most relevant ele- 
ments observable in the X-ray spectral range between 1.8 and 
10 keV, namely silicon, iron and nickel. 

Our analysis shows that, within the cool-cores of bright 
nearby clusters, metal abundances of Si, Fe and in a few in- 
stances even Ni can be measured to a high precision. Indeed the 
precision is so high, particularly for Fe, that we need to intro- 
duce a cross-calibration uncertainty of 3% to reconcile measure- 
ments secured with different EPIC experiments. Thanks to the 
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Table 6. Supernovae type la fraction, / = S NIa/(S NIa + SNcc), obtained for different theoretical SNe models from the observed 
Si/Fe and Ni/Fe abundance ratios. SNIa and SNcc models are the same as in Table|5] 



mekal model 


SNIa 


SNcc 


W7 


W70 


WDD1 


WDD2 


WDD3 


CDD1 


CDD2 


No06,S 


0.117-0.136 


0.110-0.244 












No06,T 


0.139-0.141 


0.132-0.272 












CL04,S 




0.149-0.193 




0.168 -0.384 


0.134-0.289 




0.154-0.348 


CL04,T 




0.182 -0.320 




0.204-0.313 


0.165 -0.327 




0.188 -0.302 


W95,S 






0.207 - 0.254 


0.129 -0.322 


0.102-0.237 


0.226 - 0.239 


0.118-0.289 


W95,T 








0.143 -0.338 


0.113-0.250 




0.131 -0.304 



apec model 


SNIa 


SNcc 


W7 W70 


WDD1 


WDD2 


WDD3 


CDD1 


CDD2 


No06,S 


0.096-0.109 




0.107 - 0.270 


0.086 - 0.200 




0.098 - 0.243 


No06,T 






0.130-0.301 


0.104-0.225 




0.119-0.272 


CL04,S 




0.221 - 0.606 


0.144-0.345 


0.174-0.262 


0.239-0.517 


0.133-0.313 


CL04J 




0.269 - 0.543 


0.179 -0.388 


0.145 -0.300 


0.290 - 0.453 


0.166-0.354 


W95,S 




0.172 -0.556 


0.197-0.286 




0.186-0.564 


0.189 -0.258 


W95,T 




0.191 -0.574 


0.147-0.301 




0.207 -0.517 


0.141 -0.272 



high statistical quality of our data we find evidence for some in- 
trinsic scatter in the element abundance with Fe around 20% and 
Si and Ni about 30%, abundance ratios are also characterized by 
roughly the same intrinsic scatter. 

The relatively modest scatter around mean values of Si, Fe 
and Ni abundances and Si/Fe and Ni/Fe ratios indicates that, 
whatever the process responsible for the enrichment of the ICM, 
it works rather similarly in all objects. Even in the somewhat ex- 
traordinary case of the Centaurus cluster with its anomalously 
large metal content, the abundances relative to iron are per- 
fectly compatible with those of the other clusters indicating the 
presence of similar enrichment processes. Uncertainties on the 
abundance estimates associated to the specific choice of spectral 
model, i.e. 2T rather than 4T appear to be negligible (i.e., within 
the statistical uncertainties for Si and Ni, and of the same order 
of the EPIC cross-calibration differences for iron). 

The estimate of the abundances from our spectra may be 
thought of, with some simplification, as a two step process: in 
the first step the equivalent width of a given line is estimated; 
in the second step the equivalent width is converted into an 
abundance assuming the temperature derived by fitting the con- 
tinuum. While the first step is relatively straight-forward and, 
at least for K-shell lines, leaves little room for ambiguity, the 
second, requiring estimates of collisional ionization and transi- 
tion probabilities is far more prone to differences. We have veri- 
fied that the two plasma codes available within XSPEC provide 
somewhat different estimates for the metal abundances, the dif- 
ferences are not hu ge and qualitatively similar to those foun d 
by other workers dSander s & Fabian 200^ Ide Plaa et al. I l2007h . 
however, given the high quality of our abundance estimates, they 
provide the dominant source of indetermination for Si (10%), Ni 
(16%) and Ni/Fe (22%) and contribute substantially in the case 
of Fe (4%) and Si/Fe (%6). 

We find that the Si/ Fe abundance ratio e stimated for cool- 
core clusters (this work; iTamura et alj|2004l) . shows a remark- 
ably good agre ement with the values found from samples of 
galaxy groups (Rasmus sen & Ponmanl 120071) and X -ray lumi- 
nous elliptical galaxies dHumphrev & Buotd2006l) . This similar- 
ity favors a common enrichment scenario for luminous elliptical 
galaxies and the cool-cores of groups and clusters, i.e. a similar 
mix of SNIa and SNcc. 



We have used our estimates of the Si/Fe and Ni/Fe abun- 
dance ratios to constrain the relative contribution of type la and 
core-collapsed supernovae to the enrichment process. To this end 
we have compiled a list of SNIa and SNcc yields from the liter- 
ature. Simply by inspecting our yields (see Tab. I A. 1 1 and IA.2I ) 
it is rather obvious that the large differences, in the tens of %, 
are bound to have a non-negligible impact on our estimates. The 
Fe normalized yields of Si and Ni both feature a scatter that is 
larger than the indetermination in the corresponding observed 
abundance ratios. Under the rather simplistic assumption of a 
20% indetermination in the yields we find that, /, the fraction 
of la to la plus core-collapsed supernovae cannot be reconciled 
with or 1, in other words we need both la and core-collapsed 
supernovae to produce the observed ratios. This result appears 
to be rather solid, an indetermination of at least 50% is required 
to allow / = and than only for one specific combination of la 
and core-collapsed supernovae, an indetermination of more than 
70% to have / = 1. 

Going beyond the qualitative statement that both SNIa and 
SNcc contribute to the enrichment of the ICM is quite hard. The 
accepted range for /, the ratio of type la SN to all SNe, con- 
volved over all the SNIa and SNcc combinations providing valid 
results, is rather large, from 10% to 40%. Similarly we can say 
that: roughly more than half (48-79%) of the Fe is produced by 
SNIa; less than half of the Si (14-49%) is produced by SNIa and 
Ni is very poorly constrained (26-88%). Increasing the indeter- 
mination on the yields beyond 20% will of course further enlarge 
the allowed ranges and result in even weaker constraints. We 
reiterate that the limiting factor is not the X-ray abundance es- 
timates, mean values can be constrained rather well, even when 
allowing for systematics associated to instrument calibration and 
thermal emission model uncertainties, but the large indetermina- 
tion in the theo retical yields . A sim ilar point was made some ten 
years ago by Gibson et alj d 1997b . a decade later the substan- 
tially improved observational constraints stand out in stark con- 
trast with the lack of any similar advancement on the theoretical 
side. 

Our analysis shows how difficult it is to determine the rel- 
ative contribution of type la and core-collapsed supernovae to 
the enrichment of the ICM from global cool-core measurements. 
An alternative approach is to measure radial profiles of metal 
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abundances. In this case the variation of an abundance ratio such 
as Si/Fe, where the two elements are produced in different pro- 
portions by SNcc and SNIa, can be interpreted as evidence for 
a variation of the contribution of one type of SN with respect 
to th e other. This kind of analysis has been attempted in the 
past, iFinoguenov et all d2000l 1200 ll) using ASCA data, found 
that for most rich clusters in his sample the Si/Fe rapidly in- 
creased when moving from small to large radii. Lat er work using 
XMM- Newton data did not confirm these findings. iTamura et alj 
(2004) find a flat Si/Fe ratio for a sample of 19 clusters. Amongst 
the few objects in comm on between the ASCA and XMM- 
Newton samples a re A3 1 1 2 lFinoguenov et alj d2000t) and A2052 
IFinoguenov et al J d2001l) for which ASCA data found Si/Fe 
radial gradients and XMM-Newton derived flat Si/Fe profile s. 
Recent analysis of SUZAKU data dSatoet all l20"07bl |2009albl: 
iKomivama etaT]|2009l) on a handful of groups and poor clusters 
leads to Si/Fe profiles that are consistent with being constant 
out to at least O.lrigo and in some instances to 0.2rigo- Since at 
0. lri8o the avera ge Fe abundance excess is ab out 1/2 of what it is 
within 0.03ri8o (Leccardi & Molendi 2008a), this seems to rule 
out the possibility that the transition from central excess to flat 
Fe profile might be associated to a change in SN type mix. 

Future work on XMM-Newton, Chandra and SUZAKU data 
may change somewhat the situation by providing new abun- 
dance ratio measurements beyond O.lriso. However, given the 
considerable di fficulties involved in estima ting abundances at 
large radii (e.g. iLeccardi & Molen di 2008a), it is by no means 
clear if and by how much our estimates will improve. A ma- 
jor advancement will be provided by the first mission carrying a 
micro-calorimeter, most likely the Japanese ASTRO-H. The ten 
fold increase in resolution at the Si line will allow to reduce the 
background by an order of magnitude for the H-like Si line and 
by a factor of a few for the He-like Si triplet. This will translate 
into measures of Si abundances out to at least 0.3rigo for a good 
number of nearby clusters. 



8. Summary 

We have performed a detailed study of the Si, Fe and Ni abun- 
dances in the cores of 26 nearby cool-core clusters. Our work 
may be divided into two main parts, the first relates to the mea- 
sure of the abundances and the second to the estimate of the 
relative contribution of SNIa and SNcc to the ICM enrichment 
process. Regarding the first part our main result may be summa- 
rized as follows. 

- We find that systematic uncertainties associated to the dif- 
ferent spectral modeling, namely 2T versus 4T model, are 
below a few per cent (2 - 3%). 

- We find evidence for a 3% systematic error between Si and 
Fe abundance measures secured with the three EPIC detec- 
tors. The Ni abundances do not have the statistical quality to 
investigate such small systematic errors. 

- We have verified that the mekal and apec plasma codes avail- 
able in XSPEC give somewhat different abundances values. 
Given the high photon statistics of our spectra they con- 
tribute significantly to the indetermination of Fe (4%) and 
Si/Fe (6%) and are the dominant source of indetermination 
from Si (~ 10%), Ni (~ 15%) and Ni/Fe (~ 20%). 

- The final Si, Fe and Ni abundance distributions as well as 
the Si/Fe and Ni/Fe abundance ratios distributions of the 
sample show only moderate spreads (from 20% to 30%) 
around their mean values (exact values are reported in Tab. 



H|. These nearly constant distributions suggest similar ICM 
enrichment processes at work in all cluster cores. 

- We find a remarkably uniformity in the observed Si/Fe abun- 
dance ratio ranging from X-ray luminous elliptical galaxies 
to the cool-cores of groups and clusters. This tell us that, 
whatever the real proportion between different SNe types 
may be, the enrichment process of the hot gas associated to 
elliptical galaxies is likely the same in isolated ellipticals, 
dominant galaxies in groups and brightest cluster galaxies in 
clusters. 

Regarding the second part, we have used our estimates of the 
Si/Fe and Ni/Fe abundance ratios to constrain the relative contri- 
bution of SNIa and SNcc to the enrichment process of the ICM. 
We have considered a suite of 6 SNcc yields and 7 SNIa yields 
and, unlike previous studies, we have included both uncertainties 
associated to the observed abundance ratios and the theoretical 
yields for Si, Fe and Ni. Under the assumption that the indeter- 
mination on the yields currently available in the literature is of 
the order of 20%, we find that: 

- the SNIa iron-mass-fraction, overall permitted range is 
0.48 - 0.79; 

- the SNIa silicon-mass-fraction and nickel-mass-fraction 
ranges are 0.14 - 0.49 and 0.26 - 0.88, respectively; 

- the number ratio of SNIa over the total number of SNe, /, 
spans from 0.10 to 0.38; 

- the dominant source of uncertainty in the estimate of £ are 
errors on theoretical yields, conversely, the choice of spectral 
code has almost no impact. 

Our analysis shows that the large uncertainties on the cur- 
rently available yield tables prevent any precise estimate of the 
relative contribution of SNIa and SNcc, and that all that can re- 
ally be said with some certainty is that they both concur to the 
ICM enrichment process in the cool-core regions. An alternative 
approach to estimate the relative contribution of SNIa and SNcc 
is to measure the difference between abundance ratios, such as 
Si/Fe, in the core and in the outer regions of individual systems. 
While early attempts have not provided definitive results, a more 
robust analysis of currently available data may furnish important 
constraints. The coming into operation of the first space-borne 
micro-calorimeter, most likely ASTRO-H, will undeniably allow 
a significant step forward. 
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Fo r core-collap sed sup ernovae we use the yields 

from iNomotoetalJ d2006l) . IChieffi & Limongil d2004) and 
IWooslev & Weaver! ( 1 19951) . In these cases we have first averaged 
the yields over the Initial Mass Function (IMF) of the stellar 
population as follows: 



d™ (l+X)dm 



(Al) 



where M,(m) is the z-th element mass p roduced by a star of 
mass m. We consider two IMF: a standard Salpeter (1955) func- 
tion with exponent x = 1 -35, and , a top-heavy function with 
x = 0.95 (I Arimoto & Yoshii|[l987l thereafter AY), which pre- 
dicts a relatively larger number of massive stars. The integration 
mass range is betwee n 10 and 50 solar masses in agreement with 
iNomotoet all (120061) . 

The derived yields of Si, Fe and Ni in units of solar masses 
per supernova explosion are reported in Table IA.1I For core- 
collapsed supernovae we report results for the standard Salpeter 
and the AY top-heavy IMF. 

Inspection of Table fA. 1 I reveals how much yields differ from 
one model to the other. Fe is characterized by the smallest scat- 
ter, 10% for SNIa and 18% for SNcc, Si is less constrained 27% 
for SNIa and 19% for SNcc, Ni features the largest variations, 
46% for SNIa and 41% for SNcc. Note also that the choice of 
IMF produces only a modest impact on the yields. 

The derived SNIa yie lds of Si and Ni relative to iron in 
lAnders & Grev esse ( 1989) sol ar unit s are reported in Table lA.2l 
Ratios in iGrev esse & Sauval dl998l) units can be easily com- 
puted by multiplying both the Si/ Fe and Ni/Fe v alues in Table 
IA.2l bv 0.675, whereas ratios in Lodders (2003) units are ob- 
tained bjnimltigl^ing Si/Te by 0.65 and Ni/Fe by 0.68, and, ra- 
tios in lAsplund et al.l d20 05) units by multiplying Si/Fe by 0.67 
and Ni/Fe by 0.72 (for th e solar Ni we have used the updated 
value bv lScott et al.ll2009l) . 

The derived yields of Si and Ni relative to iron in solar units 
for SNcc for different progenitor metallicities and IMFs are re- 
ported in Table IA.2I Conversion factors for different solar units 
systems are the same as given above. 

The yields relative to Fe for all o ur SNcc models ar e 
also listed in a recent review paper (Wer ner et alj [2008). 
Unfortunately some of the values reported in that paper are in- 
cor rect, in some cases, by a la rge amount, i.e. the Ni yields for 
the IWooslev & Weaver] d 19951) models. We have corresponded 
with the authors who are now considering publishing a revised 
version of their tables. 



Appendix A: 

We have computed yields for Si, Fe and Ni for different models 
of la and core-collapsed supernovae. 

The ones for SNIa originate from two ph ysically distinct sets 
of models taken from llwamoto et alJ dl999h . namely the slow 
deflagration (W7 and W70) and the delayed detonation (WDDs 
and CDDs) explosion models. We have computed silicon, iron 
and nickel yields for a given SNIa model by summing up over 
all stable isotopes the synthesi zed amount of Si, Fe and Ni per 
SN event given in TablelA~T1of Iwamot o et alJ (fl999l) . 
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Table A.l. Yields of elements in solar masses for different SNe models. Type la y ields were calculated for each element for the W7, 
W70, WDD1, WDD2, WDD3 and C DD1, CDD2 models inllwamoto et atl(ll999l) . SN core-collapsed yields (including types lb and 
Ic) were calculated from the results in lNomoto et al.l (l2006f) integrating over a progenitor initial mass function. We have considered 
two different initial mass functions, the Salpeter function with index 1.35 and the AY top-heavy function with index 0.95, and we 
have integrated between 10 and 50 solar masses. The initial metallicity of the progenitors is Z=0.02 (i.e. solar). 



SNIa 


Element 


W7 


W70 


WDD1 


WDD2 


WDD3 


CDD1 


CDD2 


Si 


0.16 


0.14 


0.27 


0.21 


0.16 


0.28 


0.20 


Fe 


0.76 


0.77 


0.67 


0.79 


0.87 


0.65 


0.83 


Ni 


0.13 


0.10 


0.04 


0.06 


0.07 


0.04 


0.06 



SNcc 

Element Nomoto etal. (2006) Chieffi&Limongi (2004) Woosley&Weaver(1995) 



Salpeter top-heavy Salpeter top-heavy Salpeter top-heavy 



Si 


0.113 


0.124 


0.163 


0.183 


0.127 


0.132 


Fe 


0.090 


0.091 


0.132 


0.135 


0.109 


0.106 


Ni 


0.005 


0.005 


0.014 


0.013 


0.014 


0.013 



Table A.2. Yields of Si and Ni relative to Fe in solar units o f lAnders & Grevesseld 19891) for the same SNe models used in Tab. IA.ll 



SNIa 


Ratio 


W7 


W70 


WDD1 


WDD2 


WDD3 


CDD1 


CDD2 


Si/Fe 
Ni/Fe 


0.55 
4.20 


0.48 
3.23 


1.07 
1.42 


0.69 
1.85 


0.48 
2.09 


1.13 
1.35 


0.63 
1.85 



SNcc 

Ratio Nomoto etal. (2006) Chieffi&Limongi (2004) Woosley&Weaver(1995) 

Salpeter top-heavy Salpeter top-heavy Salpeter top-heavy 
Si/Fe 3.28 3.58 3.22 3.56 3.07 3.26 

Ni/Fe 1.43 1.38 2.71 2.43 3.24 2.99 



